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Abstract 

Beauty production in events containing two muons in the final state has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
114 pb~^. A low transverse- momentum threshold for muon identification, in com- 
bination with the large rapidity coverage of the ZEUS muon system, gives access 
to almost the full phase space for beauty production. The total cross section for 
beauty production in ep collisions at \/i = 318 GeV has been measured to be 
c"tot(ep — >■ hbX) = 13.9 ± 1.5(stat.)l4;3(syst.) nb. Differential cross sections and a 
measurement of bb correlations are also obtained, and compared to other beauty 
cross-section measurements, Monte Carlo models and next-to-leading-order QCD 
predictions. 
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1 Introduction 



The production of beauty quarks in ep collisions at HERA provides a stringent test of 
perturbative Quantum Cliromodynamics (QCD), since the large 6-quark mass {rrib ~ 
5 GeV) gives a hard scale that should ensure reliable predictions in all regions of phase 
space, including the kinematic threshold. Especially in this region, with 6-quark transverse 
momenta comparable to or less than the 6-quark mass, next-to-leading-order (NLO) QCD 
calculations in which the (massive) b quarks are generated dynamically are expected to 
provide accurate predictions [1,2]. 

The cross section for beauty production has been measured in pp collisions at the SppS [3] 
and Tevatron colliders [4], in 77 interactions at LEP [5,6], and in fixed-target ttN [7] 
and pN [8] experiments. Most results, including recent results from the Tevatron, are in 
good agreement with QCD predictions. Large discrepancies are observed in some [5] of 
the results from 77 interactions at LEP. 

In most of the previous measurements of beauty production at HERA, beauty events 
were selected by requiring the presence of one or more jets, tagged by a muon or electron 
from the semi-leptonic decay of one of the b quarks [9-13], or by tracks originating from 
the secondary decay vertex of beauty hadrons [14]. This restricts the measurements to b 
quarks with high transverse momentum (pr)- 

This paper reports measurements of beauty production via the reaction ep — bbX — >■ 
yU/iX' using the ZEUS detector at HERA. The dimuon final state yields a data sample 
enriched in bb pairs, and with strongly suppressed backgrounds from other processes. This 
allows low muon-pT (p^) thresholds to be applied without any jet requirements, and gives 
access to a larger region of phase space, especially towards lower transverse momenta of 
the b quarks. 

Conceptually, the analysis is similar to the HI and ZEUS analyses of beauty in D*fi final 
states [15,16], with three significant differences. The larger branching ratio yields higher 
statistics, so that differential cross sections can be measured. The wider rapidity coverage 
and very low pt threshold allow the extraction of the total beauty cross section with little 
extrapolation. The low charm background in the dimuon final state, partially due to the 
harder b fragmentation, allows measurements of bb correlations, testing the influence of 
higher-order contributions on the perturbative calculations. 



2 Experimental set-up 

The data sample used in this analysis corresponds to an integrated luminosity C = 
114.1 ± 2.3 pb"\ collected with the ZEUS detector from 1996 to 2000. In 1996-97, 
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HERA provided collisions between an electroiu beam of Ef, = 27.5 GeV and a pro- 
ton beam of Ep = 820 GeV, corresponding to a centre-of-mass energy ^/s = 300 GeV 
(£300 = 38.0 ± 0.6 pb~^). In 1998-2000, the proton-beam energy was Ep = 920 GeV, 
corresponding to ^/s = 318 GeV (£318 = 76.1 ± 1.7 pb~^). 

A detailed description of the ZEUS detector can be found elsewhere [17]. A brief outline 
of the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [18], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angl^ 
region 15° < ^ < 164°. The transverse-momentum resolution for full-length tracks was 
(^{Pt)/pt = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [19] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/ \^ for electrons and 
a{E)/E = 0.35//E for hadrons, with E in GeV. 

The muon system consisted of rear, barrel (R/BMUON) [20] and forward (FMUON) [17] 
tracking detectors. The B/RMUON consisted of limited-streamer (LS) tube chambers 
placed behind the BCAL (RCAL), inside and outside the magnetised iron yoke surround- 
ing the CAL. The barrel and rear muon chambers covered polar angles from 34° to 135° 
and from 135° to 171°, respectively. The FMUON consisted of six planes of LS tubes 
and four planes of drift chambers covering the angular region from 5° to 32°. The muon 
system exploited the magnetic field of the iron yoke and, in the forward direction, of two 
iron toroids magnetised to 1.6 T to provide an independent measurement of the muon 
momentum. 

Muons were also detected by the sampling Backing Calorimeter (BAG) [21]. This detector 
consisted of 5200 proportional drift chambers which were typically 5 m long and had a wire 
spacing of 1 cm. The chambers were inserted into the iron yoke of the ZEUS detector 
(barrel and two end caps) covering the CAL. The BAG was equipped with analogue 



^ Electrons and positrons are both referred to as electrons in this paper. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as 77 = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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readout for energy measurement and digital readout for muon tracking. The former was 
based on 2000 towers (50 x 50 cm^), providing an energy resolution of ~ 100%/ V^. The 
digital information from the wires allowed the reconstruction of muon trajectories in two 
dimensions {XY in barrel and YZ in end caps) with an accuracy of a few mm. 

3 Principle of the measurement 

Events with at least two muons in the final state were selected. Two principal event classes 
contribute to the beauty signal to be measured. The first consists of events in which the 
two muons originate from the same parent b quarljfl, e.g. through the sequential decay 
chain b — > cfiX — > sfi^X'. These yield unlike-sign muon pairs produced in the same 
event hemisphere and with dimuon invariant masses of m'^J^^ < 4 GeV (i.e. a partially 
reconstructed 5-meson mass). The second class consists of events in which the two muons 
originate from different beauty quarks of a bb pair. These can yield both like- and unlike- 
sign dimuon combinations, depending on whether the muon originates from the decay of 
the primary beauty quark, or from a secondary charm quark. In addition, B^B^ mixing 
can dilute these charge correlations. Muons from different b quarks will predominantly 
be produced in different hemispheres, and tend to have a large dimuon mass. 

An important background contribution arises from primary charm-quark pair production 
where both charm quarks decay into a muon. This yields unlike-sign muon pairs only, with 
the two muons produced predominantly in opposite hemispheres. Since this background 
is too small to be measured directly from the dimuon data, it was normalised to the charm 
contribution as determined from the ZEUS D* + fi sample [16] which has a similar event 
topology and covers a similar though somewhat more restricted kinematic range. 

Other backgrounds yielding unlike-sign muon pairs include heavy quarkonium decays and 
Bethe-Heitler (BH) processes. In contrast to muons from semileptonic decays, muons from 
these sources are not directly accompanied by hadronic activity, thus giving an isolated 
muon signature. 

Beauty production is the only source of genuine like-sign muon pairs. Background contri- 
butions to both like- and unlike-sign combinations include events in which either one or 
both muons are false, i.e. originate from — > or tt — > /x decays in flight or are misidenti- 
fied hadrons. Studies [22] have shown that the charges of such false-muon pairs are almost 
uncorrelated, i.e. the contributions to the like- and unlike-sign dimuon distributions are 
almost equal. The difference between the unlike- (iV^ta) like-sign {N\^^^) distribu- 
tions is thus essentially free from false-muon background, without the need to simulate 

^ Unless stated otherwise, throughout this paper, the term b quark includes b. 
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this background using Monte Carlo (MC) methods. Once the background contributions 
from open charm (iVcharm), J /''P and other heavy vector mesons (AVm) and Bethe-Heitler 
(A^bh) are known, this difference can be used to measure the beauty contribution N^i^^^ 
according to the formula 

Ar,5^,, = (A^l,, - A^La - (A^charm + iVvM + iVBH)) X _ f (1) 

where the last term refers to the unlike-sign (A^^) and like-sign (A^^^) beauty contributions 
predicted by the MC. Small corrections to this procedure will be explained in Section [71 
The beauty signal is hence extracted from the difference between the unlike- and like-sign 
samples. 

The like-sign false-muon background can then be obtained from the data by subtracting 
the MC like-sign beauty contribution, properly scaled to the measurement, from the total 
like-sign sample, while the unlike-sign background is a simple reflection of the like-sign 
background. This method to obtain the false-muon background contributions will be 
referred to as the subtraction method. 

Since one of the goals is the determination of the total beauty production cross section in 
ep collisions, events from deep inelastic scattering (DIS), where the photon virtuality, Q^, 
is larger than 1 GeV^, and photopro duct ion (Q^ < 1 GeV^) were not explicitly separated. 

The average cross sections obtained from the two different running periods (y^ = 300 
and 318 GeV) are all expressed in terms of a single cross section at =318 GeV. This 
involves a typical correction of +2%. 



4 Event selection and reconstruction 
4.1 Trigger selection 

The data were selected online by means of a three-level trigger system [17,23] through an 
inclusive "or" of four different trigger channels: 

• a muon reaching the inner B/RMUON chambers and matched to a minimum ionis- 
ing energy deposit (MIP) in the CAL or any muon reaching the outer B/RMUON 
chambers (muon channel); 

• a reconstructed D meson candidate {D* channel [24], plus similar chains for other 
charm mesons [25]); 

• two jets (dijet channel [11]); 



4 



• a scattered-electron candidate in the CAL (DIS channel [12]). 

For part of the data taking, the requirements on the DIS and dijet channels were loosened 
in the presence of any muon in the inner B/RMUON chambers. The non-muon trig- 
gers were used to gain geometric acceptance for regions not covered by the B/RMUON 
chambers, and to evaluate the efficiency of the muon triggers. Owing to this redundancy, 
the trigger efficiency for dimuon events with reconstructed muons from beauty was high, 
80 ± 4%. 

4.2 Event selection 

The large mass of a bb pair, at least ~ 10 GeV, usually leads to a significant amount of 
energy deposited in the more central parts of the detector. To suppress backgrounds from 
false- muon events and charm, a hadronic transverse energy cut 

Et>8 GeV 

was applied, where 

{E!lr'^^° no scattered electron 

^e^>io° _ j^e^ with scattered electron. 

The transverse energy was calculated as Ef^^^^° = Ej 0.>io°(-E'i sin^j), where the sum runs 
over all energy deposits in the CAL with the polar angle above 10°. The latter restriction 
is imposed to remove proton-remnant effects. If detected, the energy of the scattered 
electron (E^) was subtracted. The detection criteria for the scattered electron were the 
same as in a previous publication [16]. 

Various tracking requirements were imposed [22], the most important of which was that 
the reconstructed longitudinal vertex position should be consistent with an ep interaction, 
l^vtxl < 50 cm. 

4.3 Muon selection 

Muons were reconstructed offline using an inclusive "or" of the following procedures: 

• a muon track was found in the inner B/RMUON chambers. A match in position and 
angle to a CTD track was required. In the bottom region, where no inner chambers 
are present, the outer chambers were used instead. If a match was found to both inner 
and outer chambers, a momentum-matching criterion was added; 



5 



• a muon track was found in the FMUON chambers. Within the CTD acceptance, a 
match in position and angle to a CTD track was required and the momentum was 
obtained from a combined fit of the CTD and FMUON information. Outside the 
CTD acceptance, candidates well measured in FMUON only and fitted to the primary 
vertex were accepted; 

• a muon track or localised energy deposit was found in the BAC, and matched to a 
CTD track, from which the muon momentum was obtained. In the forward region 
of the detector, a MIP in the calorimeter was required in addition in order to re- 
duce background related to the proton beam or to the punch-through of high energy 
hadrons. 

Most muons are within the geometric acceptance of more than one of these algorithms. 
The overall efficiency is about 80% for high-momentum muons (more than 2-5 GeV, 
depending on 77). 

Two different kinematic selections were made. In the barrel region, the requirement that 
the muons reach at least the inner muon chambers implies a muon transverse momentum 
(j9^) of about 1.5 GeV or more. In order to have uniform kinematic acceptance, a cut 

> 1.5 GeV 

was therefore applied to all muons (selection A). 

In the forward and rear regions, lower px muons can be detected, although with somewhat 
higher background. To cover the largest possible phase space for the intended measure- 
ment of a total beauty-production cross section, the px cut was lowered to 

> 0.75 GeV 

for high-quality muons [22], i.e. muons seen by more than one muon detector and/or 
confirmed by a MIP in the CAL (selection B). For other muons satisfying all previously 
listed criteria, the cut > 1.5 GeV was retained to keep the background low. Selection 
A is thus a subset of selection B. 

At least two such muon candidates were required per event. No explicit cut on the muon 
angle was applied for either selection. The angular coverage of the muon chambers, BAC 
and CTD gives continuous useable acceptance in the pseudorapidity region 

-2.2 < r/^ < 2.5 . 

To suppress events with ambiguous matches between CTD tracks and muon chamber 
segments as well as genuine dimuons from prompt light-meson decays (e.g. p — /i/i), a 
dimuon invariant mass (m^^') cut of 

m^^ > 1.5 GeV 
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was applied. This implies a minimum opening angle between the two muons. 

Events with a very forward and a very backward muon candidate, a topology not favoured 
for the beauty signal, were removed by a cut on the difference in pseudorapidity of the 
two muon candidates of 



Muon candidates with badly measured momentum (predominantly from false-muon back- 
grounds) were suppressed using the imbalance between the transverse momenta of the 
muons 



An additional cut with a similar scope as the initial Ex cut was applied on the fraction 
of the total transverse energy carried by the muon pair 



The reason for the distinction of the two different dimuon mass regions will be explained in 
Section [71 This iJ^-fraction cut removes events where the hadronic activity is, respectively, 
very high (false-muon background) or low (quarkonia and Bethe-Heitler). 

Cosmic-ray muons were removed by discarding events with back-to-back muon candidates 
and events in which the average calorimeter timing differs by more than 10 ns from the 
nominal collision time. Large cosmic showers were removed using the BAG total energy 
and number of BAG muon segments. 

A sample of 4146 dimuon events was obtained using selection B. Selection A retained 
about two thirds of these events. 

4.4 Muon isolation 

Muons from semileptonic decays are usually not isolated, i.e. they are normally accompa- 
nied by hadrons originating from the fragmentation and decay of the parent heavy quark 
and from other hadronic activity in the event. Hadronic activity in the detector was re- 
constructed using a combination of both track and calorimeter information [26] referred 
to as energy-flow objects (EFOs). The difference in azimuth angle and pseudorapity, A0 
and At], was calculated between each EFO and each muon candidate in the event. The 
total transverse energy, Ji 2, deposited in a cone of AR = y^A^^TAr/^ < 1 around each 
muon flight direction was calculated by summing over all relevant EFOs, excluding the 
other muon. Since usually either both (beauty signal and open charm) or neither (elas- 
tic J/ip, Bethe-Heitler, etc.) of the muons arise from semileptonic decays, the quadratic 



1]' 



-r/'^^l < 3.0. 



iW-pfP\)/ifT+P'T)<0.7. 




for m''^ < 4 GeV 



for m/^^ > 4 GeV. 
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sum = If + 12 of the two energy sums was found to yield the best sensitivity to 
distinguish between the two cases. 

5 Background and event simulation 

In order to measure the beauty signal, several background contributions to the selected 
data sample were evaluated: 

• the background from open charm decays not originating from beauty; 

• the background from quarkonium states not originating from open beauty {J/ip, 
T, ...), produced in elastic or inelastic collisions; 

• the background from Bethe-Heitler muon pair production; 

• the background from false muons. 

Monte Carlo simulations of beauty and charm production were performed using the gen- 
erators Pythia [27] (for events with < 1 GeV^) and Rapgap [28] (for > 1 GeV^). 
These simulations include the direct photon-gluon fusion process {jg QQ, Q = b,c), 
flavour excitation in the resolved photon and proton (e.g. Qg — > Qg, — > Qg), and 
hadron-like resolved photon processes (e.g. gg QQ). Gluon splitting into heavy flavours 
{g QQ) in the initial or final states of light-quark events was not included in the sim- 
ulations; this contribution is, however, expected to be small [29]. 

Inelastic quarkonium production was simulated using Herwig [30], while elastic quarko- 
nia and Bethe-Heitler processes were produced using several generators including Grape 
[31]. 

The ZEUS detector response, including the transformation of MC truth level quantities 
into reconstructed quantities, was simulated in detail using a programme based on Geant 
3.13 [32]. The detector simulation for beauty and charm events includes the simulation 
of both real and false muons. 

Fake muons can be produced by hadron showers leaking from the back of the calorimeter 
or by charged hadrons traversing the entire calorimeter without interaction. In addition, 
low-momentum muons can originate from in-fiight decays of pions and kaons. Tracks 
reconstructed in the central tracker may also be erroneously associated to a signal from a 
real muon in the muon chambers. A study [29] based on pions from decays, protons 
from A decays, and kaons from and D* decays, showed that the detector simulation 
reproduced these backgrounds reasonably well. They will be collectively referred to as 
false muons. 
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Backgrounds from false muons in events not containing charm or beauty were not sim- 
ulated. They were estimated from the data using the subtraction method described in 
Section [3l 

Since the muon range in dense material (effective momentum threshold) and the muon 
detector efficiencies were imperfectly simulated, corrections to the MC were determined 
[22] using an independent data set consisting of isolated J/ip and Bethe-Heitler events. 
Tabulated as a function of and r/'*, these corrections were applied to MC events on an 
event-by-event basis. 

6 Theoretical predictions and uncertainties 

The MC programs described earlier, based on leading-order (LO) matrix elements with 
the addition of parton showers (PS) to obtain higher-order topologies, were used for the 
acceptance corrections. These programs are expected to describe the shapes of differential 
distributions, but not necessarily their normalisation. For quantitative comparisons with 
QCD, next-to-leading-order (NLO) predictions are used. 

QCD calculations in which h quarks are treated as massless particles [33] are not applicable 
in the kinematic range relevant here. Calculations based on CCFM parton-evolution 
schemes [34] , also called kx factorisation, do not yet exist with full NLO implementation. 
Fixed-order NLO calculations with massive h quarks were therefore chosen as the reference 
predictions. 

The NLO FMNR program [1] evaluates parton-level cross sections for beauty in 7p colli- 
sions (photoproduction) in the fixed-order massive mode, for both pointlike and hadron- 
like photon couplings to the heavy quarks. The Weizsacker- Williams (WW) approxi- 
mation with an effective Q^ax cutoff of 25 GeV^ (~ ml) [35] was used to evaluate and 
include the DIS contribution to the cross sections, which is approximately 15%. This is 
in agreement with the DIS prediction from HVQDIS described below. 

The parton-density functions used were CTEQ5M [36] for the proton, and GRV-G-HO [37] 
for the photon. The renormalisation and factorisation scales /i were chosen to be equal 
and parametrised by /io = ^/p^^^^rn^/2, where pt is the average transverse momentum of 
the two emerging h quarks, and m^, = 4.75 GeV is the 6-quark mass. Such a scale choice is 
equivalent to the choice /xo = -E't/2 or /io = \/ + used in many jet measurements 
at the Tevatron [38] and at HERA [39], and is expected to compensate somewhat for 
uncalculated higher-order contributions [40]. An estimate of the theoretical uncertainty 
was obtained by simultaneously varying 4.5 < < 5.0 GeV and fio/2 < /i < 2/io such 
that the uncertainty was maximised. Typical uncertainties resulting from this procedure 
(e.g. for the bb total cross section) are +60% and —30%. Variations of the parton densities 
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and the strong coupling parameter, Aqcd? led to uncertainties which were much smaller 
than the uncertainties related to mass and scale variations. They were therefore neglected. 

Predictions for visible /i/i final states were obtained by linking the FMNR parton-level 
predictions to the fragmentation and decay chain provided by Pythia using the FMNR®- 
Pythia interface [41]. Additional parton showering was not appliecj^. The branching 
ratios were corrected to correspond to those obtained from the Particle Data Group 
(PDG) [43], as listed in Table [H All other parameters, including those for fragmentation, 
and the procedure to obtain their uncertainties, were the same as in an earlier analysis [16], 
and described elsewhere [41]. 

The DIS part of the inclusive cross section is also calculated using the NLO predictions 
from HVQDIS [2]. Only point-like contributions are included in this prediction. The 
parton density function used was CTEQ5F4 [36]. The renormalisation and factorisation 
scales /i were chosen to be equal and parametrised by yUo = \/Vt + + Q?' 1'^- The 
mass and scales were varied as for FMNR. A scheme for the calculation of visible cross 
sections for correlated final states, corresponding to the FMNR®Pythia interface de- 
scribed above, was not available. Therefore, DIS cross-section comparisons are limited to 
parton level, and the DIS contribution to the inclusive cross sections is included in the 
FMNR®Pythia predictions via the WW approximation. 

7 Signal extraction 

Dimuon mass and charge separation 

As motivated in Section[3], events were separated by the muon charges into like- and unlike- 
sign dimuon samples. To differentiate between muon pairs from the cascade decay of the 
same b quark and those from different b quarks, the distributions were further separated 
depending on the dimuon invariant mass: low-mass dimuons with m^^ < 4 GeV, enriched 
in muons from the same b quark, and high-mass dimuons with m'^'^ > 4 GeV, containing 
dimuons originating from the decay of different b quarks only. The dominant signal and 
background contributions to the four subsamples are summarised in Table |2j 

The resulting dimuon mass distributions for the low- and high-mass, like- and unlike-sign 
subsamples for selection B are shown in Fig. [H The MC distributions were in each case 
normalised to the data according to the procedure described in the following subsections. 
The high-mass region is already strongly beauty enriched, while the low-mass region 

* The MC@NLO approach [42], which allows the combination of NLO matrix elements with parton 
showers, is not yet available for ep interactions. 
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exhibits a significant contribution from J/ip production not originating from B liadro: 
decays. Such dimuon pairs tend to be isolated. 



Dimuon isolation cuts 

To reduce this J /if) contribution, as well as corresponding contributions from , T and 
Bethe-Heitler processes, a non-isolation requirement was applied, based on the fact that 
muons from semileptonic decays are accompanied by hadrons from the same decay. The 
dimuon isolation variable I^^, defined in Section [4.41 was required to exceed 250 MeV, 
safely above the noise level of the CAL. This reduces the elastic quarkonium and Bethe- 
Heitler contributions to an almost negligible level. 

Inelastic quarkonium and Bethe-Heitler events might pass the above cut because hadrons 
from e.g. the proton remnant can accidentally end up in the isolation cone. For events in 
the J/ip and ip' mass peaks, where this background is largest, the cut was therefore raised 



are called non-isolated. This additional requirement is satisfied by 3500 events from 
selection B. The other events form a complementary isolated background sample. 

Figure [2] shows the muon pt and rj distributions for non- isolated unlike-sign dimuon pairs, 
combining the low- and high-mass samples. The remaining contribution from J/ip, Bethe- 
Heitler, etc. processes was normalised to the isolated background sample. The charm 
contribution is small and was normalised to the charm signal in the D* + fi sample [16] 
as outlined in Section [21 The different contributions to Fig. [2] are listed in Table [31 

Signal evaluation 

The beauty signal and false-muon background were obtained using the procedure de- 
scribed in Section [31 However, some further corrections are needed. Events from unlike- 
sign background sources, such as charm, which have been reconstructed as like-sign 
dimuon events due to false muons, are included both in the false-muon background esti- 
mation and in the MC samples. To avoid double counting, this (very small) contribution 
is subtracted from the MC samples. False muons in the beauty MC are considered as 
part of the signal. 

^ The term B hadron includes b baryons. 



to 2 GeV. 



In summary, dimuons fulfilling the relation 




2.0 GeV for m^'^' E [2.9, 3.25] GeV or m^'^' e [3.6, 3. 
0.25 GeV otherwise 



75] GeV 
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The signal extraction procedure according to Eq. ([T]) relies on the unlike- and like-sign 
false-muon background contributions being equal. A dedicated false-muon background 
study [22] revealed a small residual excess of unlike-sign over like-sign background. This 
excess was corrected for using a multiplicative correction factor, Ocorr, of 1.02 for the high- 
mass and 1.06 for the low-mass dimuon sample. The beauty fraction was thus determined 
using a modified version of Eq. ([T]) 

/ iV^- + N^- \ 

Nbl^,, = (A^data - «corr-iVLa " (A^charm + Nyu + ATbh)) X -—^ ^ . (2) 

A total of 1783 of the 3500 non-isolated events from selection B were found to originate 
from beauty, corresponding to a beauty fraction of 51%. 



8 Systematic uncertainties 

The main sources of systematic uncertainty for the measurement of visible cross sections 
are described in this section, in approximate order of importance. The numbers in paren- 
theses refer to the specific case of the inclusive visible cross section of Section [9l Bin-by-bin 
uncertainties were evaluated for the differential distributions where possible and appro- 
priate. They are mostly similar to those derived for the inclusive visible cross section. 
Additional uncertainties introduced by the extrapolation to quark-level cross sections are 
discussed in Section |9l 

• Muon efficiency correction. The muon efficiency, including the efficiency of the 
muon chambers and of the MUON-CTD matching, is known to about 7% from a study 
based on an independent muon sample, and from the variance of the cross section when 
information from different muon detectors is used independently [22]. Conservatively, 
it is assumed to be fully correlated between the two muons (±15%). 

• Normalisation of charm background. The transfer of the normalisation of the 
charm contribution from the D*fi analysis [16] to this analysis involves the following 
uncertainties: the statistical error of the fit of the charm contribution, ±10%; the 
inclusive branching ratio c ^ fi, ±10%; the acceptance uncertainty due to charm 
fragmentation and decay spectra, ±10%; the fragmentation fraction c D*^, ±6%; 
the branching ratio D*^ Kittt, ±3%; and the use of all muon detectors (this 
analysis) versus the use of the barrel and rear muon chambers only [16], ±10%. The 
influence of the correlation between the fitted beauty and charm fractions in the D*fi 
analysis [16] was found to be negligible. The normalisation of the charm contribution 
was varied by 21% according to the resulting combined uncertainty (±12%). 
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Normalisation of the Bethe-Heitler, J/t}), etc. backgrounds. The normalisa- 
tion of the residual non-isolated contributions from Bethe-Heitler, charmonium, and 
T production was varied by ±50% (±10%). 

False-muon background. As a cross check for the determination of the false-muon 
background by the subtraction method, the probability of a reconstructed hadron to be 
misidentified as a muon was obtained from an inclusive dijet MC sample and tabulated 
as a function of pt and rj. Starting from a data sample with selection cuts identical to 
the present analysis, except that only one muon candidate was required, false-dimuon 
events were created by assuming a suitable additional hadron to be identified as a 
muon according to this tabulated probability. After corrections for trigger efficiency, 
and for the contribution from one false and one true muon obtained directly from the 
h and c MC, an independent background prediction was obtained [44]. It agreed very 
well in both normalisation and shape with that of the default subtraction method, thus 
confirming the method. Since the uncertainty on this background is already implicitly 
contained in the statistical error of the subtraction method, no explicit additional 
uncertainty was assigned. 

h spectral shape uncertainty and hb correlations. It was checked that the h- 
quark spectra from Pythia and Rapgap agree well with the corresponding spectra 
from the NLO predictions described below [29]. To estimate the effect of variations of 
this shape, and of effects of variations of the hb correlations for different topologies on 
the efficiency, the efficiency was evaluated using the Pythia direct contribution only, 
or doubling the non-direct contributions (±4%/— 12%). 

B^B^ oscillations. The B^B^ oscillation parameter was varied by 8%. This includes 
the uncertainties of the mixing implementation in the MC models used (±4%). 

Other b MC model uncertainties. This includes the uncertainty of the procedure 
used to account for differences of the branching ratios in the signal MC and Table [H 
the uncertainty from b fragmentation, and from the shape of the lepton spectrum from 
b decays (± 10%). 

c spectral shape uncertainty and cc correlations. The direct and non-direct 
fractions for the charm background MC were varied in the same way as for beauty. 
The effect on the signal was small (±0%/— 4%). 

Trigger efficiency. The uncertainty on the trigger efficiency was estimated by com- 
paring the efficiencies for muon and non-muon triggers in data and MC (±5%). 

Other uncertainties. Other uncertainties include the variation of the like-/unlike- 
sign ratio for the false-muon background by 3% (±3%), the variation of the isolation 
cuts by up to 500 MeV (±2%), the variation of the Et cut (energy scale) by 3% 
(±2%), the variation of the cuts (magnetic field uncertainty) by 0.3% (< 1%). 
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The total systematic uncertainty (+25%/— 28%) was obtained by adding the above con- 
tributions in quadrature. The uncertainties related to the background normalisation and 
the b and c spectral shape uncertainties were applied at a bin-by-bin level where relevant, 
while the others were added globally. A 2% overall normalisation uncertainty associated 
with the luminosity measurement was not included. 

9 Total bb cross section 

As a first step towards the extraction of the total cross section for bb production, a visible 
cross section was extracted for the maximum possible region in muon phase space allowed 
by the preselection and the detector acceptance (selection B). The criterion that the muon 
detection probability should be at least about 30% per muon leads to the following phase 
space definition at truth level: 

• —2.2 < 1] < 2.5 for both muons; 

• Pt > 1-5 GeV for one of the two muons; 

• Px > 0.75 GeV for the other muon, as well as p > 1.8 GeV for rj < 0.6, or {p > 2.5 
GeV or > 1-5 GeV) for rj > 0.6. 

This cross-section definition refers to only one pair of muons per event. If there are more 
than two muons, muons directly originating from B hadron decays are taken preferentially 
to form the pair. A visible cross section for dimuon production from beauty decays in 
this phase space 

avis(ep ^ bbX fifiX') = 55 ± 7(stat.)+J^(syst.) pb (3) 

was obtained. This cross section includes muons from direct 5-hadron decays, and indirect 
decays via intermediate charm hadrons or r leptons. The two muons can either originate 
from the same b quark, or from different quarks of the bb pair. Muonic decays of kaons, 
pions or other light hadrons were not included. 

The measured cross section is larger than, but compatible with, the FMNR®Pythia 
NLO prediction 

t^vis,NLo(ep ebbX efifiX) = 33l^^(NLO)+3(frag. © br.) pb, (4) 

where the first error refers to the uncertainties of the FMNR parton-level calculation, and 
the second error refers to the uncertainties related to fragmentation and decay. 

The visible cross section was then translated into the total cross section for beauty produc- 
tion. The effective branching fraction of a bb pair into at least two muons is 6.3% [27,43]. 
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The probability (acceptance) for such a muon pair to be in the kinematic range of the 
measured visible cross section, evaluated from the beauty MC sample, is about 6% on 
average. Defining as the maximum of the two 6-quark transverse momenta after 
parton showering, and |Cfe|™° as the minimum of the modulus of the rapidity (not pseudo- 
rapidity) of the two quarks, this probability ranges from 3% for p™^^ = GeV to 9% all§ 
PtT ~ 10 GeV, for iCbl™™ < 2. The acceptance is almost independent of rapidity within 
this rapidity range, which covers 90% of the total bb phase space. It drops sharply at 
larger rapidities. Thus, only 10% of the total beauty contribution in the region |Cb|™™ > 2 
remains unmeasured. The small dependence of the acceptance on the transverse momenta 
of the b quarks is due to the low muon-momentum threshold, in combination with the 
large 6-quark mass and the three-body decay kinematics. Sensitivity down to = GeV 
is obtained. 

In summary, the combined probability for a bb pair to yield a muon pair in the visible 
kinematic range (6.3%x6%=0.38% on average) is quite small, but varies by less than a 
factor 3 over 90% of the total phase space. Furthermore, it is almost entirely determined 
by quantities measured [43] with good precision at e'^e" colliders. These include the 
branching fractions listed in Table [H the 6-fragmentation functions, and the B hadron — > 
fiX decay spectra. It was checked that all of these are well reproduced by the MC after the 
application of branching-ratio corrections. The 6-quark pt and rapidity spectra predicted 
by the Pythia and Rapgap generators were found to agree with those from FMNR 
and HVQDIS to within 15% [29]. Furthermore, the quasi-uniformity of the acceptance 
explained above implies that the dependence on details of the simulation of the bb topology 
is rather weak. The MC can therefore safely be used for the extraction of the total cross 
section for beauty production. 

The normalisation of the Pythia + Rapgap MC prediction for the beauty contribution 
had to be scaled up by a factor 1.84 to agree with the dimuon data. Applying this 
measured scale factor to the total Pythia and Rapgap cross sections, the total cross 
section for bb pair production in ep collisions at HERA for -^i = 318 GeV was determined 
to be 

fTtot(ej5 ^ bbX) = 13.9 ± 1.5(stat.)+t3(syst.) nb, (5) 

where the first uncertainty is statistical and the second systematic. In addition to the 
uncertainties described in Section [HI this includes an error of 5% from the uncertainties 
of the spectral shape mentioned above, and an error of 6% from the variation of the 
branching ratios, added in quadrature. 

The total cross section predicted by next-to-leading-order QCD calculations was obtained 
in the massive approach by adding the predictions from FMNR [1] and HVQDIS [2] for 

^ At even larger the acceptance rises further, but the fraction of events is smah. 
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less than or larger than 1 GeV^, respectively. The resulting cross section for s/s = 318 
GeV 

al^^{ep-.bbX) = 7.5ttinh 

is a factor 1.8 lower than the measured value, although compatible within the large uncer- 
tainties. The corresponding cross section from FMNR only using the Weizsacker- Williams 
approximation to estimate the DIS contribution is 

CW(ep ^ bbX) = 7.8tti nb, (6) 
in agreement with the more exact FMNR+HVQDIS calculation. 

The fact that the comparisons between data and theory yield the same ratio at the visible 
level (Eqs. ©/(iD): 

pdata/NLO _ y+0.7 
^vis ~ ' -1.1 ! 

and at quark level (Eqs. ([5])/ ([6])): 

jjdata/NLO q+0.8 

confirms the validity of the extrapolation procedure used. 

Figure [3] shows a comparison of the measured total cross section to cross sections and 
theoretical predictions from the D* + /i final state obtained by ZEUS in earlier measure- 
ments [16]. Although not fully inclusive, these measurements are closest in phase space 
to the measurement presented here. Qualitatively, they show the same trend of the cross 
sections being higher than, but consistent with, the corresponding QCD predictions. The 
somewhat larger deviations reported in similar D* + measurements by HI [15] are not 
supported. 



10 Differential cross sections and bb correlations 

Selection A was used for the measurement of visible differential cross sections because 
a uniform kinematic acceptance is more relevant than maximal phase-space coverage. 
Correspondingly, at truth level, the phase space was restricted to: 

• > 1.5 GeV for both muons 

• -2.2 <r]'' < 2.5 . 

The backgrounds were again normalised as described in Section [71 The signal-extraction 
procedure was the same as for the inclusive visible cross section, except for being applied 
bin by bin. Bin-dependent systematic uncertainties were calculated wherever possible. 
The resulting cross sections for the differential pj. and t]^ spectra are shown in Figs. H] 
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and O Very good agreement is observed with the Pythia+Rapgap predictions scaled 
by the same factor 1.84 that was measured for the total cross section. Apart from the 
normalisation, the leading-order plus parton-shower (LO+PS) approach yields a good 
description of the corresponding physics processes within the entire accessible phase space. 
This confirms the applicability of these MC models for acceptance calculations. 

A comparison of the measured cross sections to the absolute FMNR®Pythia NLO QCD 
predictions is also shown in Figs. H] and [51 Again, good agreement in shape is observed, 
with a tendency to underestimate the data normalisation consistent with the observations 
from the total cross section. A potential trend for increasing data/theory deviations 
towards low px and/or high 77, suggested by some previous measurements [10, 12], is not 
supported. 

To provide a more detailed look at the correlations between the two b quarks, the re- 
constructed dimuon mass range was restrictecill to m'^'^ > 3.25 GeV. This reduced the 
contribution of dimuons from the same quark to an almost negligible level. The corre- 
sponding data distribution for A(j) between the two muons is shown in Fig. [61 Figure [7] 
shows the resulting differential cross section, where the mass cut was replaced by the re- 
quirement that the two muons originate from different b quarks. The distribution is well 
described by the FMNR®Pythia NLO QCD predictions within the large uncertainties 
resulting from the subtraction method (Eq. ([2])). 

11 Hadron- and parton- level cross-sections 

In order to compare to previous ZEUS results using other final states [11, 12], expressed 
in terms of parton-level cross section^ differential in p^, similar cross sections were also 
extracted. 

The first step was the extraction of visible cross sections for B hadrons in different pt 
ranges. For this purpose, the data sample used for the measurement of the total beauty 
cross section (selection B) was split into two subsamples, with m'^'^ > 3.25 GeV and m'^'^ < 
3.25 GeV. As motivated in the previous section, the m'^'^ > 3.25 sample is dominated 
by muons from different b quarks, with correlations between the two quarks which are 
reasonably understood. Thus, two measured B hadrons are present in each event. To 
estimate their transverse momenta, the quantitity 

E^"" =p^^ + I^ 

^ While the mass separation value of 4 GeV described earlier was optimised such that all dimuons from 
the same h quark contribute to the low-mass sample, including dimuons from b ^ ipf decays, the value 
3.25 GcV was chosen to optimise the separation power for dimuons from same and different b quarks. 

^ HI results have not been published in this form. 
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is evaluated for each muon, where I'^ is the cone transverse energy described in Section 
14.41 This variable is found to be strongly correlated to the parent B hadron transverse 
momentum at high p^, where the additional energy from fe-quark fragmentation to the B 
hadron compensates the loss due to the unreconstructed neutrino from the semileptonic 
decay. At Pr ^ ^b, this correlation is diluted by the effect of the 5-hadron mass and the 
corresponding decay kinematics. Figure [H](a) shows the expected 5-hadron pt spectra for 
three bins in E^'', < E^'"" < 5 GeV, 5 < E^'^ < 10 GeV, and 10 < E^'^ < 40 GeV. 
Reasonably distinct i?-hadron px regions are probed. The corresponding visible cross 
sections are shown in Fig. [9](a). 

A similar procedure was applied to the subsample with m^^ < 3.25 GeV. In this sample, 
the muons originate mainly from the same b quark, therefore only one B hadron has been 
measured. Due to branching ratios and decay kinematics, the cross section is smaller, but 
the absence of like-sign muon pairs from the same b quark leads to a smaller uncertainty 
from the subtraction method. Therefore, the precision of the measurement is comparable 
to that from the high-mass region. Furthermore, the subtraction method reduces the 
influence of the residual contribution of muons from different b quarks. Thus, the measured 
cross sections are almost completely insensitive to bb correlations. 

The E^^^ variable is redefined to 

—Pt+ -'high 

where is the transverse momentum of the dimuon system added vectorially, and J^jgj^ 
is the isolation of the higher p^ muon only, to avoid double counting. The correlations 
to the B hadron pt are similar to the high- mass case (Fig. [Uj^b)), enabling them to be 
combined later on. The resulting visible 5-hadron cross sections are shown in Fig. [9](b). 

For both subsamples, agreement is found with the FMNR®Pythia predictions, consistent 
with the conclusions obtained earlier. 

The second step is to extrapolate these cross sections to 6-quark level. For comparison 
with previous measurements, the cross sections were restricted to photoproduction. Each 
of the 5-hadron visible cross sections is translated into a differential cross section 4nr in 
the pseudorapidity range |?7t| < 2 [11] with photon virtuality Q"^ < 1 GeV^ and inelasticity 
0.2 < y < 0.8, using the FMNR(8>Pythia predictions. Each cross section is quoted at the 
mean value for events satisfying the cuts for the corresponding E^^^ bin. The results 
are shown in Fig. [TOl 

The cross sections derived from the low- and high-mass subsamples (same and different 
b quarks) are in agreement, and were combined to give a single cross section for each 
p^ value. The maximum possible correlation of the systematic errors is assumed for this 
combination. 
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The resulting combined cross sections are compared to theory and previous measurements 
in Fig.[TTl They are consistent with these previous measurements, and extend the mea- 
sured range to lower p\^. Predictions at NLO [1] and predictions from a LO /cy-factorisation 
approach [34] yield an equally good description of the data. 

12 Conclusions 

The total cross section for beauty production in ep collisions at ^/s = 318 GeV has 
been measured for the first time using an analysis technique based on the detection of 
two muons, mainly from semileptonic beauty decay. The almost complete phase-space 
coverage combined with the weak dependence on details of the bb event topology allowed 
a reliable extraction of the total beauty production cross section, with acceptance down to 
= GeV, and a direct comparison to NLO QCD predictions. The predictions are lower 
than the observed cross sections, but compatible within the uncertainties. Differential 
cross sections in pi^, rj'^, and A0^^ were also measured. Shapes predicted by Monte Carlo 
models incorporating leading-order matrix elements followed by parton showers agree 
well with the data. NLO QCD predictions agree in shape with both the data and the 
LO-I-PS predictions, but are again somewhat lower than the data, in agreement with the 
observation from the total cross section. The angular correlations between final-state 
muons from different b quarks, reflecting the correlations between these parent quarks, 
are described by the NLO QCD predictions. Measurements of cross sections for muon 
pairs from the same or from different B hadrons yield similar and compatible results. A 
comparison with previous measurements through the extrapolation to differential cross 
sections at 6-quark level shows reasonable agreement, and extends these measurements 
down to lower p^. 
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channel 


effective branching fraction w/o B^B^ mixing 


b ji" direct 


10.95 ± 0.27 % 


6 ^ /i"^ indirect 


8.27 ±0.40 % 


b ^ fi~ indirect 


2.21 ±0.50 % 


all b /i± 


21.43 ±0.70 % 


bb (diff. bs) 


2.42 ±0.17 % 


bb fi^fi^ (diff. 6s) 


2.18 ±0.14 % 


b fJ'^fJ'" all 


2.40 ± 0.16 % 



Table 1: Effective branching fractions used for cross-section determinations. The 
indirect contributions include cascade decays into muons via charm, anticharm, 
and J/ip. The additional effect of B^B^ mixing (x = 0.1283 ± 0.0076^ is not 
included. 





unlike-sign ±/± 


like-sign ± + / 


low inv. mass 
m^^ < 4 GeV 


muons from same b, 

muons from J/'0, ip', 
and false-muon background 


false-muon background, 

and small contribution of 
muons from different b 


high inv. mass 
m^f, > 4 GeV 


muons from different b, 

muons from cc, T, BH, 
and false-muon background 


muons from different b 

and false-muon background 



Table 2: Classification of events using dimuon mass and charge correlations. 
The main contributions to each class are listed; the most relevant is indicated in 
bold face. 
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process 


muon candidates 


beauty 


2382 


charm 


629 


quarkonia and BH 


281 


false muon 


1281 


data 


4574 



Table 3: Number of muon candidates contributing to Fig. 2: unlike-sign non- 
isolated dimuons. 4^74 muons correspond to 2287 events. 
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Figure 1: Dimuon mass distributions of unlike- sign dimuon pairs from selection 
B (see text) in the (a) low-mass and (c) high-mass subsamples, as well as like-sign 
dimuon pairs in the (b) low-mass and (d) high-mass subsamples. The same vertical 
scale has been chosen for the like- and unlike-sign subsamples, with different bin 
sizes for the high- and low-mass regions. The expected contributions from different 
processes are also shown. The false-muon background was obtained from the data 
using the subtraction method described in Section{^ Due to this method, the total 
prediction for like-sign pairs agrees with the data by definition. 
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Figure 2: (a) Muon transverse momentum and (b) muon pseudorapidity distri- 
bution from both high- and low-mass dimuon pairs in the non-isolated unlike-sign 
sample. Two muons are entered for each event. The expected contributions from 
different processes are also shown. Due to the subtraction method, the statistical 
error of the prediction for the false-muon background is comparable in absolute size 
to that of the data. 
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Figure 3: Comparison of measured cross sections to NLO QCD predictions. The 
bb cross section from this analysis (top) is compared to both measured and predicted 
borb cross sections obtained in the ZEUS D*fi analysis [16] for the photoproduction 
regime (middle line) and DIS (lower line). The NLO calculations in the D*fi 
analysis used a slightly different set of parameters. Using the parameters detailed 
in SectionlBi the central value of the photoproduction cross-section prediction would 
increase by about 20%. 
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Figure 4: Cross-section da/dpf^ for muons from b decays in dimuon events with 
pt^ > 1.5 GeV and —2.2 < rj^ < 2.5 for both muons. Two muons contribute for 
each event. The data ( solid dots ) are compared to the scaled sum of the predictions 
by the LO+PS generators Pythia and Rapgap (histogram) and to the NLO QCD 
predictions from FMNR®Pythia (band). 



ZEUS 





45 


a. 




a. 


40 


?^ 

■o 


35 


"b 




■o 


30 




25 




20 




15 




10 




5 








ZEUS 114 pb-1 
(PYTHIA+RAPGAP) x 1.84 
NLO QCD (FMNRSPYTHIA) 



-1.5 -1 -0.5 0.5 



1.5 2 2.5 



Figure 5: Cross-section da/drj^ for muons from b decays in dimuon events with 
pf^ > 1.5 CeV and —2.2 < t]^ < 2.5 for both muons. Two muons contribute for 
each event. The data (solid dots) are compared to the scaled sum of the predictions 
by the LO+PS generators Pythia and Rapgap (histogram) and to the NLO QCD 
predictions from FMNR®Pythia (band). 
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Figure 6: Distribution of the azimuthal distance A(f) between the two muons 
in dimuon events with > 1.5 GeV for both muons, and m^^ > 3.25 GeV. The 
expected contributions from different processes are also shown. Due to the subtrac- 
tion method, the statistical error of the prediction for the false muon background is 
comparable in absolute size to that of the data. 
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Figure 7: Cross-section da/dA(j)^'^ for bb events in which the muons origi- 
nate from different b quarks, with p^ > 1.5 GeV and —2.2 < rj^ < 2.5 for both 
muons. The data (solid dots) are compared to the scaled sum of the predictions by 
the LO+PS generators Pythia and Rapgap (histogram) and to the NLO QGD 
predictions from FMNR®Pythia (band). 
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Figure 8: Distribution of the true pr of the parent B hadron for muons from (a) 
different b quarks or from (b ) the same b quark, for the three E^'^ bins indicated in 
the figures. For the definition of see Section [771 
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Figure 9: Visible cross section for parent B hadrons from events containing two 
muons satisfying the cuts for the total cross-section measurement, and in which 
both muons originate from a different (a) or from the same b{b) quark (b), in three 
bins of p^. There are two entries per event for (a), and one entry per event for 
(b). The data (solid dots) are compared to the scaled sum of the predictions by 
the LO+PS generators Pythia and Rapgap (histogram) and to the NLO QCD 
predictions from FMNR®Pythia (band). 



31 



o 

i 5- 
OS 

Q 3h 



ZEUS 



T 



ZEUS 114 pb"" 
o different b 

* same b 

• average 
NLO QCD 



2 - 



10 



I 

15 20 

(GeV) 



Figure 10: Data/NLO ratio for cross sections from different b quarks (open 
circles ) compared to measurements from the same h quark ( stars ) and their average 
(filled circles). The value for each E^'^ (or p^) interval (0-5,5-10,10-40 GeV) is 
quoted at the median Pt of the parent b quarks in events satisfying all detector 
level cuts (4-7,8.0,14-0 GeV)- The three points for each p\. value are shown slightly 
shifted in p\ for clarity. 
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Figure 11: Differential cross section da/dp^ of this analysis (stars) compared 
to previous ZEUS measurements (other symbols), FMNR NLO QCD predictions 
(band), and predictions from the factorisation approach (thick line). 
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